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In QCD, the similarity of multiplicity distributions is violated i) by the running of 
the strong coupling constant a s and ii) by the self-similar nature of parton cascades. 
, It will be shown that the data collapsing behavior of P n onto a unique scaling 

^ ■ curve can be restored by performing the original scaling prescription (translation 

and dilatation) in the multiplicity moments' rank. 

m 

1 Similitude 

The notion of scaling is hardly new. One of the earliest scaling arguments 
, dates back to 1638 when Galileo Galilei published his., infamous masterpiece 

entitled "Dialogues Concerning Two New Sciences ".u Among other funda- 
mental observations he examined the principle of similitude, the elementary 
properties of similar physical/biological structures. Galileo realized that the 
strength S of a bone increases in direct proportion to its cross-sectional area 
(S ~ I 2 , ii I is the linear size), whereas the weight of a bone increases in 
direct proportion to its volume (W ~ I 3 ). Thus, there will be a characteristic 
point where a bone has insufficient strength to support its own weight: the 
intersection point of the quadratic and cubic curves denoting the strength and 
weight of a bone, respectively. This general engineering consideration implies 
that terrestrial bodies can not exceed a certain maximum size. The classical 
scaling argument of Galileo teaches us an important lesson: the physical laws 
are not invariant under a uniform change of the size of macroscopic objects. 
The gravitational force, governed by Newton's constant Gjv with dimension 



o 
o 



X 



of (mass) -2 , inevitably leads to the breakdown of dilatation symmetry. 

Classical scaling principles of the above sort are based on the key assump- 
tion that the physical bodies or processes are uniform, filling an interval in 
a smooth, continuous fashion. In the example given by Galileo, the strength 
of a bone was assumed to be uniformly distributed over the cross-sectional 
area with its weight having a similar uniformity. This is a major limitation of 
the principle of similitude because such assumptions are not necessarily accu- 
rate. In reality a vast number of biological and physical systems, the so-called 
fractals, exhibit highly irregular appearance as the result of their self-similar 
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structure. Let us consider a well-known example, the architecture of the hu- 
man lung. If one unfolds the cca. 300 million air sacs of the self-similar 
bronchial tree and merges them into one continuous flat surface, its area will 
be as large as a tennis court. This anomalous surface-to- volume ratio can 
not be explained by classical scaling arguments based on Galileo's principle of 
similitude. Only the modern, more powerful scaling ideas of fractal geometry 
can properly characterize self-similar geometric forms. 

2 Similar Distributions 

Is it meaningful to speak about the concept of similarity with regard to multi- 
plicity fluctuations? Of course, yes. Counting the number of particles created 
in a certain collision process, one of the most basic observables is the distri- 
bution of counts: the multiplicity distribution P n . It is a discrete distribution 
but at high energies we can approximate P n by a continuous probability den- 
sity f{x) either via P n rs f(x — n) or via P n rs f^~n +1 f( x ) where f(x) is 
called similar if it satisfies 



with A > being a scale parameters In multiparticlc physics one often sets 
c = 0, A = (n(s)) and uses P n (s) ~ f(x = n, s) to approximate the shape 
of P n (s) measured at different collision energies s. Then Eq. (1) means that 
expressing the multiplicities n in units of (n(s)), the properly rescaled data 
points, preserving normalization, fall onto the universal curve ip{ z ) which 
depends only on the dimensionless ratio z = n/ (n(s)). This behavior is called 
KNO scaling after the work of Koba, Nielsen and Olesen.B Two years earlier 
it was obtained by Polyakov,Q too. Sometimes people try to improve on the 
scaling via shifting the multiplicity distributions by a factor c(s) ~ 1. Usually 
this number is interpreted as the average of produced leading particles. 

Can we extend the similarity property (1) for multiplicity distributions 
P n (5) measured at different bin-sizes 5 in phase space? Not quite. The exper- 
imental data collected in the past 15 years or so revealed a dominant feature 
of multiplicity fluctuations: in a wide range of collision energies, bin-sizes, and 
for a large variety of reaction types, the observed fluctuation pattern proved 
to be self- similar. This so-called intermittent behavior manifests through the 
power-law dependence of thx^njormalized factorial moments of P n (S) as the 
resolution scale S is varied J3i3'l0 whereas Eq. (1) expresses the constancy of 
normalized moments. The breakdown of the similarity feature Eq. (1) due 
to self-similar multiplicity fluctuations is analogous to the incompatibility of 
Galileo's principle of similitude and the properties of fractal geometric forms. 
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3 Scaling and Quantum Mechanics 



As we have seen previously, an obvious reason of the breakdown of dilatation 
symmetry of physical laws is the appearance of explicit scales, such as the 
masses of macroscopic bodies or of elementary particles. But there is another 
source of non-scale-invariance, related to the properties of the quantum me- 
chanical vacuum. In quantum mechanics the physical vacuum is a polarizable 
medium. Virtual pairs of charges are always present as quantum mechanical 
fluctuations whose effect can not be switched off. They partially screen or 
antiscreen a test charge. Therefore its effective value depends on the distance 
or energy scale at which it is measured. In other words, the effective coupling 
strength is running in quantum theory. This fundamental effect has important 
consequences for multiplicity fluctuations, too: the various scaling behaviors 
inevitably break down at certain energy and resolution scales. For example, 
in c + c~ annihilation the .s-depcndencc of the QCD coupling constant a a can 
not be compensated by a suitable change of A and c in Eq. (1). Thejjnniiing 
of a s is expected to cause violation of KNO scaling at high energiesJjlJElj 

4 New Multiplicity Scaling Law 

The multiplicity moments (n q ) provide another very useful representation of 
the information encoded in P n . Our variable in this case is the rank q. Is it 
meaningful to perform a scaling transformation of type (1) in the moments' 
rank? If so, what kind of dynamics yield a shifting or rescaling in q-space? 
The Mellin transform of a probability density f(x) is defined by A4{f(x); q} = 
J Q x q ~ x f(x) dx and it provides the moment (for simplicity we make 

use of P„ ~ f(x = n)). Via the functional relation 

^M(f(x);^\=ML r f(x^;q^ (2) 

one can introduce translation and dilatation in the moments' rank q by per- 
forming the transformation f(x) — ► x r f[x^j of the probability density f(x) 
approximating the shape of P n . The above scaling relation in .M-space is our 
main concern in the remaining sections. 

5 Dilatation in Mellin Space 

The most important source of dilatation in Mellin space is related to QCD. 
In higher-order pQCD calculations, allowing more precise account of energy 
conservation in the course of multiple parton splittings, the natural variable of 
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the multiplicity moments is the rescaled rank qj instead of rank q itself.EHIll 
Here j(a s ) is the so-called QCD multiplicity anomalous dimension. Because 
of the running of the strong coupling constant a s , it is inevitable to adjust 
an energy dependent scale factor in Mellin space if we want to arrive at data 
collapsing of P n {s) onto a universal scaling curve. 
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Figure 1. a) The MLLA pQCD prediction Eq. (3) for scaling violation: V( 2 ) with fi = 5/3 
(dots) asymptotically evolves to V"( z ) with fj, = 1 (solid line), b) The scaling behavior is 
recovered in the logarithmic scaled multiplicity x = ln(Dz) according to Eq. (4). 

Let us consider in detail the shape change of P n (s) in c + c~ annihilation.! 
Taking into account MLLA corrections responsible for energy-momentum 
conservation in parton jets, the analytic form of the KNO scaling function 
becomes a gamma distribution in the power-transformed variable z M : 

4>{z) = JW^" 1 exp ( - [DzY) (3) 

where k = 3/2, Af = fiD^ k /L(fc), fi = (1 - 7)- 1 w 5/3 and D is a scale 
parameter depending on 7(a s ); 7 « 0.4 at LEP-1 energy. Thus, the MLLA 
calculation predicts violation of KNO scalingjj see Fig. la, since /i varies with 
collision energy s due to the running of a s . Note, however, that data collapsing 
can be restored in a simple manner using logarithmic scaling variable; for the 
KNO function Eq. (3) we get 

ip(x) = fie^(kfxx-e^ x )/T(k), x = \n(Dz). (4) 

Because only the exponent /i and scale parameter D of (3) are expected to 
depend on collision energy s through the variation of 7(a s ), data collapsing is 
recovered by plotting fi~ 1 ip(fix) as displayed in Fig. lb. The scale change in 
logarithmic multiplicity is governed by the multiplicity anomalous dimension 
of QCD, which sets the scale in Mellin space, too - see oarx basic relation (2). 
This type of scaling of P n (s) is called log-KNO scaling^ since one observes 
the behavior of type (1) but now the distribution of logarithmic multiplicity 
turns out to be similar. 
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In c + e~ annihilation the breakdown of ordinary KNO scaling at high 
energies is only expected to arise. In hh collision, however, this proved to be a 
dominant feature of observations already in the mid-80s when the exploration 
of SPS energies started. With the log-KNO law in our hands it is challenging 
to test its validity using real data. The violation of Eq. (1) is most visible for 
multiplicities measured by the E735 Collaboration.!^ The full phase space 
multiplicity distributions were obtained in pp and pp collisions at cm. ener- 
gies y/s = 300, 546, 1000 and 1800 GeV at the Tevatron collider. 

At Tevatron energies, bimodal 
shapes of the distributions show up 
having shoulder structure - like at 
SPS. It was argued liij that the low 
multiplicity regimes are influenced 0.01 
mainly by single parton collisions > 
and exhibit KNO scaling, whereas 
the large- n tails of the distributions 
are influenced more heavily by dou- 0001 
ble parton interactions and violate 

L -10 12 

(1) considerably. This part of the 4 

data sets was analyzed in log-KNO Fi § ure 2 - Log-KNO scaling of the E735 data, 
fashion and, as shown in Fig. 2, scaling holds with good accuracy. Our (still 
preliminary) investigation suggests that double parton collisions yield a scale 
change not only in multiplicity but in the multiplicity moments' rank as well, 
whereas single parton collisions do not produce the latter effect. 




6 Translation in Mellin Space 

The other major source oLlhfi breakdown of Eq. (1) is the self-similarity 
of multiplicity fluctuations .EfO'Q This can be observed through the power-law 
scaling C q cx 5~ iPq of the normalized moments C q = (n q ) / (n) q of P n (S) as 
the bin-size S in phase space is varied (we neglect the influence of low count 
rates). The simplest possibility is the monofractal fluctuation pattern. Then, 
the so-called intermittency exponents tp q are given by ip q — ^2(9 — 1) and the 
anomalous fractal dimensions D q — cp q /(q — 1) are g-independent, D q = D^. 
The normalized moments C q of P n (S) take the form 

C q =A q [C 2 }"- 1 for q>2, (5) 

with coefficients A q independent of bin-size S. Eq. (1) is obviously violated 
since one measures (5-dependent second moment, C2 oc 6~ D2 . 
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In the restoration of the similarity feature (1) for self-similar fluctuations, 
the basic idea is the investigation of the higher-order moment distributions 



P n ,r = n r P n /(n r ). Their moments are 



,9+r 



)/(r 



i.e. the moments 



of the original P n are transformed out up to r-th order by performing a shift 
in Mellin space, see Eq. (2). For r = 1, the normalized moments of the first 
moment distribution P n ^i are found to be C 9 ,i = C q +\/[C2\ q in terms of the 
original C q and comparison to Eq. (5) yields C Qt i = A q+ i for monofractal mul- 
tiplicity fluctuations. Since the coefficients A q are independent of bin-size <5, 
we see that monofractality yields not only power-law scaling of the normal- 
ized moments of P n but also data collapsing behavior of the first moment 
distributions P n .\ measured at different resolution scales S. The effect of low 
multiplicities (Poisson noise) can be taken into account via the study of facto- 
rial moment distributions P n>r = rJ r l P n / '(n' r ') and their factorial moments. 

Increasing the rank of the moment 
distributions allows the restoration of 
data collapsing behavior in the pres- 
ence of an increasing degree of multL 
fractality of self-similar fluctuations O 
This feature is best seen for random 
multiplicative cascades which interpo- 
late between monofractaJs and fully 
developed multifractals.ta The fluctu- 
ations give rise to the log-Levy law 
having a characteristic parameter, the 
Levy index < v < 2. The moments 
obey the same structure as in Eq. (5) 
with exponent (q v - q)j(T - 2). For 
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Figure 3. D2 >r /D 2 for a few Levy index v. 

v = this gives back the monofractal case, whereas the upper limit of the Levy 
index, v = 2, corresponds to the log-normal law resulting from fully developed 
multifractal fluctuations. Log- normal distributions exhibit two remarkable 
properties: the higher-order moment distributions are also log-normals (form 
invariance), further, they differ from each other only up to a change of scale 
(scale invariance). Hence, for fully developed multifractals it is impossible to 
arrive at data collapsing behavior via translation in Mellin space, no matter 
how large is r, because the normalized moments remain unaltered. In the 
other limit, monofractals produce data collapsing already for r = 1. Fig. 3 
illustrates the changing degree of fractality with increasing r through the 
variation of the ratio Di^l^i '■ the larger is the value of the Levy index v, 
the harder is to arrive at P>2. r = 0. The fixed-point at v — 2 is apparent (the 
mathematically disallowed values v > 2 bring farther away from scaling) . 
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7 Summary 

In QCD, the similarity feature Eq. (1) of multiplicity distributions P n breaks 
down. For P n (s), the running of the strong coupling constant a s gives rise to 
the scale breaking. For P n (S), the self-similar nature of multiplicity fluctua- 
tions in parton jets results in the violation of KNO scaling. (Due to running 
coupling effects, self-similarity itself also breaks down at very small 8). But 
if we switch from P n to (n q ) 1 it turns out that both QCD effects can be com- 
pensated by a suitably chosen shifting and rescaling in the moments' rank q. 
That is, in order to arrive at data collapsing of the multiplicity distributions 
onto a unique scaling curve, the original similarity prescription (translation 
and dilatation) is still satisfactory, only the mathematical representation of 
fluctuations should be changed from distributions to their moments - in the 
intermittency era this is the dominant practice, anyway. The functional rela- 
tion Eq. (2) tells everything about how the scaling behavior manifests for the 
distributions themselves: [x q ^) corresponds to /(x M ) and therefore log-KNO 
scaling of the form pT 1 /(/ilnx) shows up, whereas (x q+r ) implies that the 
moment distributions x r f(x)/ (x r ) exhibit similarity. 
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pelijk Onderzoek (NWO) and the Hungarian Science Foundation (OTKA) 
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